Quantum critical spin liquids and superconductivity in the cuprates 
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We describe a new kind of quantum critical point in the context of quantum anti-ferromagnetism 
in 2d that can be understood as a quantum critical spin liquid. Based on the comparison of exponents 
with previous numerical work, we argue it describes a transition from an anti-ferromagnetic Neel 
ordered state to a VBS-like state. We argue further that the symplectic fermions capture the proper 
degrees of freedom in the zero temperature phase that is the parent to the superconducting phase in 
the cuprates. We then show that our model reproduces some features found recently in experiments 
and also in the Hubbard model. 



INTRODUCTION 

Over the last few years there has been much inter- 
est in finding new quantum critical points in the context 
of quantum anti-ferromagnetism in 2d. One motivation 
are their possible applications to the anti-ferromagnetic 
phase of the Hubbard model and to superconductivity in 
the cuprates. Strong arguments were given by Senthil 
et. al. that there should exist a quantum critical point 
that separates a Neel order phase from a valence-bond 
solid like phase In the model these authors consid- 
ered, it appears difficult to exhibit this critical point per- 
turbatively. In this paper we consider a different, simpler 
model which contains a critical point that can easily be 
studied with a perturbative renormalization group analy- 
sis. Our critical exponents agree very favorably with the 
numerical simulation of the model in [l[ carried out by 
Motrunich and VishwanathQ, and we interpret this as 
strong evidence that our model is in the same universality 
class. In the second part of the paper we take some initial 
steps toward applying this theory to superconductivity. 

For the remainder of this Introduction, wc summarize 
the motivations given in [3:| for our model. It is well- 
known that the continuum limit of the Heisenberg anti- 
ferromagnet constructed over a Neel ordered state leads 
to a non-linear sigma model for a 3-component field n(x) 
satisfying n 2 — 1. (For a detailed account of this in Id 
and 2d, with additional references to the original works 
seejij].) In Id a topological term Sg arises directly in 
the map to the continuum and affects the low-energy 
(infra-red (IR)) limit: half- integer spin chains are gap- 
less whereas the integer ones are gapped. This is the 
well-known Haldane conjecture^. It is known from the 
exact Bethe-ansatz solution of the spin 1/2 chainQ that 
the low-lying excitations are spin 1/2 particles referred 
to a spinons[7(. 

In 2d one still obtains a non-linear sigma model. The 
topological term Sg also arises but, unlike in Id, it is 
a renormalization group (RG) irrelevant operator so we 



discard it. The non-linear constraint 



= 1 renders 



captured by the following euclidean space action: 

Swf = J d 3 x Q<V* ' d n H + A(n-n) 2 ^ (1) 

dl = y; 3 ,_i dl . The fixed point is in the Wilson-Fisher 
univerality class 10] . The fixed point generalizes to an M 
component vector n and we will refer to the fixed point 
theory as off where D = d + 1. 

Senthil et. al.[l[ have given numerous arguments 
suggesting that anti-ferromagnets can have more exotic 
quantum critical points that are not in the Wilson-Fisher 
universality class. They are expected to describe for 
instance transitions between a Neel ordered state and 
a valence-bond-solid (VBS)-like phase, and some evi- 
dence for such a transition was found by Motrunich and 
VishwanathQ. See also [ll|. A large part of the litera- 
ture devoted to the "deconfined" quantum critical points 
represent the n field as 



= X f <?X 



(2) 



where a are the Pauli matrices and \ — (xi; X2) = {Xi} 
is a two component complex bosonic spinor. The con- 
straint n 2 = 1 then follows from the constraint x^X = !• 
Coupling x to a U(l) gauge field with the covariant 
derivative = <9 M — i , then by eliminating the non- 
dynamical gauge field using it's equations of motion, one 
can show that the following actions are equivalent: 



d D x dfj,n ■ dpi 



(3) 



the model non-renormalizable in 2d. However there is 
a quantum critical point in the spin system 0, 0] that is 



Senthil et. al. considered a model where x was a bo- 
son, and added an F 2 V term which makes the gauge field 
dynamical. 

The central idea of this paper is that the spinon x is 
a fermion. Numerous arguments were given in0]. First 
of all, the equivalence (|3|) is valid whether x is a boson 
or fermion. Secondly, suppose the theory is asymptot- 
ically free in the ultra-violet, which it is. Then in this 
conformally invariant limit, one would hope that the de- 
scription in terms of n or x have the same numbers of 
degrees of freedom. One way to count these degrees of 
freedom is to compute the free energy density at finite 



2 



temperature. For a single species of massless particle, 
the free energy density in 2c? is 



T = -c 3 



C(3) T 3 
2tt 



(4) 



where C3 = 1 for a boson and 3/4 for a fermion. 
is Riemann's zeta function.) (In Id the analog of the 
above is T — —cnT 2 /6, where c is the Virasoro central 



charge [ 12l |13|.) Therefore one sees that the 3 bosonic 
degrees of freedom of an n field has the same C3 as an 
N — 2 component x field. One way to possibly under- 
stand the change of statistics from bosonic to fermionic 
is by simply adding a Chern-Simons term to ([3|)[14l|. 

In analogy with the bosonic non-linear sigma model, 
since the constraint x^X — 1 again renders the model 
non-renormalizable, we relax this constraint and consider 
the action: 



J d D x (2d M x%X + 16^ 2 A \xh\ 



(5) 



where now x is an A-component complex field, some- 
times referred to as a symplectic fermion, x^X — 
J2iLiXiXi- This model may at first appear peculiar, in 
that the field x nas a Klein-Gordon action but is quan- 
tized as a fermion. However, we remind the reader that 
there is no spin-statistics theorem in 2c?. Note there is no 
gauge field ( "emergent photon" ) in the model. 

As shown in [3], the (x^x) 2 interactions drive the the- 
ory to a new infrared stable fixed point, we refer to as 
. For N = 2 in 3D the exponents were computed 
to bed 

n = 3/4, v = 4/5, P = 7/10, 6 = 17/7 (5j4 3) ) (6) 

(The definition of these exponents is given in the next 
section.) These agree very favorably with the critical ex- 
ponents found in @: v = .8 ± 0.1, j3/v = .85 ± 0.05, 
certainly within error bars. The shift down to 3/4 from 
the classical value rj = 1 is entirely due to the fermionic 

(3) 

nature of the x fields. We thus conjecture that the Sp 2 
model describes a deconfined quantum critical spin liq- 
uid. 

In the next section we summararize the results of the 
critical theory found in [3j. In section III we apply this 
model to high T c superconductors and describe agree- 
ment with some recent experimental results (l5j. 



THE CRITICAL THEORY 

The 1-loop beta function for the N component model 
in D space-time dimensions is 



dX 



(4 - D)\ + (N- 4)A 2 



(7) 



where increasing the length e corresponds to the flow 
toward low energies. The above beta function has a zero 
at 



A* — 



4- D 



N 



(8) 



Note that A* changes sign at N = 4. It is not necessarily 
a problem to have a fixed point at negative A since the 
particles are fermionic: the energy is not unbounded from 
below because of the Fermi sea. Near A* one has that 
dX/d£ ~ (D — 4) (A — A*) which implies the fixed point is 
IR stable regardless of the sign of A*, so long as D < 4. 
Arguments were given in Q that at JV" = 4, two of the x 
components reconfine into a spin field ft, though we will 
not need this here. 



Definition of the exponents for the n field 

Though the spinons x are deconfined, it is still physi- 
cally meaningful to define exponents in terms of the orig- 
inal order parameter n, which is represented by eq. 
These exponents are especially useful if one approaches 
the fixed point from within an anti-ferromagnetic phase. 
We then define the exponent rj as the one characterizing 
the spin-spin correlation function: 



(n(x)-n(O)) 



1 



| x p-a-H, 



(9) 



For the other exponents we need a measure of the de- 
parture from the critical point; these are the parameters 
that are tuned to the critical point in simulations and 
experiments: 



S x + d D x (m 2 x f X + B ■ ft) 



(10) 



Above, m is a mass and B the magnetic field. The corre- 
lation length exponent is, and magnetization exponents 
/3, 5 are then defined by 



(11) 



Above (ft) is the one-point function of the field n(x) and 
is independent of x by the assumed translation invari- 
ance. 

The above exponents are related to the anomalous di- 
mension of the field x an d the operator x^X- This leads 
to the following relations amoung the exponents: 



/3 = v(D-2 + r))/2, 6 



D+{2-rj) 
D - (2 - rf) 



(12) 



The lowest order contributions to the anomalous dimen- 
sions of the operator x*X a rise at 1-loop, and for x a t 
two loops. The calculation in [3J gives in 3D: 



2(4 - N) 
7-N '■ 



2N 2 - 17N + 33 



N 2 



11N + 28 



(13) 



3 



For N = 2 one obtains the results quoted in the intro- 
duction. 

model has the 



It was conjectured in[3j] that the Sp_ N 

' ,(3) 



same fixed point as the O^' model, so that Sp_\ is the 
3D Ising model. The exponents are in very good agree- 
ment with known Ising exponents. A shorter version of 
these results was described in 16 1 . 



SUPERCONDUCTIVITY BASED ON 
SYMPLECTIC FERMIONS 

In this section we explain how our quantum critical 
spin liquid could be relevant to the understanding of su- 
perconductivity in the cuprates, which is believed to be 
a 2 + 1 dimensional problem [17]. To do this, one must 
turn to the language of the Hubbard model. In the anti- 
ferromagnetic phase of the Hubbard model, the spin field 
n — c^ac, where c are the physical electrons. Therefore 
in applying our model to the Hubbard model, the sym- 
plectic fermion x is a descendant of the electron, so it 
can carry electric charge. Consider the zero temperature 
phase diagram of the cuprates as a function of the density 
of holes. At low density there is an anti-ferromagnetic 
phase. Suppose that the first quantum critical point is 
a transition from a Neel ordered to a VBS-like phase 
and is well described by our symplectic fermion model 
at N = 2. Compelling evidence for a VBS like phase 
has recently been seen by Davis' groupflij]; and it in fact 
resembles more a "VBS spin glass". The superconduct- 
ing phase actually originates from this VBS-like phase. 
It is then possible that the 2-component x fields cap- 
ture the correct degrees of freedom for the description 
of this VBS-like phase. These fermionic spin 1/2 spinon 
quasi-particles acquire a gap away from the critical point, 
which is described by the mass term in eq. (fT0|) . Note 
that away from the quantum critical point, the particles 
already have a gap m because of the relativistic nature 
of the symplectic fermion (20j. 

Superconductivity based on the symplectic fermion has 
some very interesting features. In the VBS-like phase the 
X-particles are charged fermions and it's possible that ad- 
ditional phonon interactions, or even the x 4 interactions 
that led to the critical theory, could lead to a pairing in- 
teraction that causes them to condense into Cooper pairs 
just as in the usual BCS theory. Recent numerical work 
on the Hubbard model suggests that the Hubbard inter- 
actions themselves can provide a pairing mechanism fljj]. 

Passing to Minkowski space, the hamiltonian of the 
symplectic fermion is 

h = J d 2 x(2a t x t a*x + 2Vx t -vx 

+ ™Vx + A (x f x) 2 ) (14) 
Expand the field in terms of creation/annihilation oper- 



ators as follows 
X(x) = 



X f (x) 



d 2 k 

d 2 k 



(«k i 

(4 



~ ik x + bu e ik xN ) 



: ik x + b] e" ik x 



(15) 



where Wk = \/k 2 + to 2 . Canonical quantization of the 
X-fields, {x^(x), <9 t x(x')} = iS(x — x')/2, leads to the 
anti-commutation relations 



(16) 



The free hamiltonian is then 

H = [ d 2 k w k (4a k + blb k ) (17) 



The minus sign in the anti-commutator of the a's 
means there are negative norm states in the free Hilbert 
space [21]. However a simple projection onto even num- 
bers of a-particles gives a unitary Hilbert space. In a po- 
tential physical realization, since the anti-ferromagnetic 
spin field n is deconfined, it is clear that the particles 
come in pairs. 

The minus sign in eq. (|16[) actually leads to a two-band 
theory. This has been seen experimentally 15] and also in 
the Hubbard model[18]. There are two kinds of spin 1/2 
particles created by a or b: ajjO) = |k) a , & k |0) = | k)t, , 
with energies e a ,b : 



H |k) a:6 = £ ,&(k) |k) a>6 
£ & (k)=w k , e a (k) = -cj k 



(18) 



Note that e a < —to and > to so there is a gap 2m. 
The density of states can easily be computed in the free 
theory. The density of states per volume is defined so 
that 



(19) 



where n is the particle number density. Using 
p^yi f d 2 k — J dee/2n, one finds 



— fb(e) for£>TO 

Z7T 



(20) 



= 



2tt 



for — TO < £ < TO 



f a (e) for £ < 



-m 



where f a ,b are temperature dependent Fermi-Dirac filling 
fractions. Interactions will tend to fill the gap. 

The last ingredient one needs is a pairing phase tran- 
sition, so let us turn to the interactions. The (x^x) 2 
interaction is very short ranged since it corresponds to 
a (5-function potential in position space. Because of the 
relativistic nature of the fields, the interaction gives rise 
to a variety of pairing interactions. There are actually 



4 



pairing interactions between the two bands. However let 
us focus on the pairing interactions within each band that 
resemble BCS pairing. If all momenta have roughly the 
same magnitude |k|, then the interaction gives the terms 
(up to factors of tt) : 



Hint — —A 



E 



-k,j 



a kj + (a — ► b) ) + 



(21) 

The overall minus sign of the interaction is due to a 
fermionic exhange statistics. Because of the overall mi- 
nus sign this is an attractive pairing interaction as in 
BCS. One difference is that in addition to the opposite 
spin pairing interactions with i ^ j, there are also equal- 
spin pairings. The quantum ground state can be further 
studied by reasonably straightforward application of the 
mean- field BCS construction [l9(. 



CONCLUSIONS 

We have shown that the 2-component relativistic sym- 
plectic fermion appears to have some of the right ingre- 
dients to explain the zero temperature phase diagram of 
the high T c cuprates. It has a quantum critical point 
that we have interpreted as a transition between an anti- 
ferromagnetic phase and VBS-like phase. Away from the 
critical point the quantum spin liquid has a 2-band struc- 
ture as in the VBS spin-glass phase [15}. It also natu- 
rally has BCS- like pairing interactions. A real test of our 
model would be a measurement of the critical properties 
of the anti-ferromagnetic to VBS spin-glass phase. The 
magnetic exponent 5 is probably the easiest to measure 
and our theory predicts 5 = 17/7. 
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